Atomic structure of Mn wires on Si(OOl) resolved by scanning tunneling microscopy 
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At submonolayer coverage, Mn forms atomic wires on the Si(OOl) surface oriented perpendicular 
to the underlying Si dimer rows. While many other elements form symmetric dimer wires at room 
temperature, we show that Mn wires have an asymmetric appearance and pin the Si dimers nearby. 
We find that an atomic configuration with a Mn trimer unit cell can explain these observations 
due to the interplay between the Si dimer buckling phase near the wire and the orientation of the 
Mn trimer. We study the resulting four wire configurations in detail using high-resolution scanning 
tunneling microscopy (STM) imaging and compare our findings with STM images simulated by 
density functional theory. 



The large magnetic moment related to a half-filled 
d-shell renders Mn atoms attractive building blocks 
for fascinating magnetic nanostructures [IH3]. Tunable 
ferromagnetism in Mn-doped semiconductors has been 
achieved for GaAs, InAs and Ge [4-7 . For silicon this 
effort has not been as successful because of strong seg- 
regation and the interstitial diffusion of Mn in the Si 
crystal during overgrowth or annealing, even though Mn 
implanted Si samples exhibit very high Curie tempera- 
tures [8]. During submonolayer deposition at room tem- 
perature (RT) however, Group III (Al, Ga, In), Group 
IV (Sn, Pb) and a few other metals (e.g. Sb and Mg) are 
known to form ID wires perpendicular to the Si dimer 
rows on the Si(100)-2xl surface [9HT7], These wires con- 
sist of metal atoms in the parallel-dimer configuration 
linking up to form atomic chains [18-20 . Recent exper- 
iments show that similar wire formation occurs for Mn 
[21]. Density functional theory (DFT) calculations pro- 
posed several possible structures for these wires [22H24] 
with one, two or three Mn atoms per Si dimer row [see 
Fig. Ha)]. 

We show that high-resolution STM imaging at RT re- 
veals a unique appearance of Mn wires. In contrast to 
other metal wires, their signature is characterized by the 
pinning of nearby Si dimers in addition to two indepen- 
dent asymmetric realizations of the Mn wire. We find 
that only an extended trimer model can explain our find- 
ings and identify a total of four wire configurations based 
on the relative orientation of the Si dimers and an asym- 
metric Mn trimer unit cell. We confirm our model using 
a sequence of STM images, where we observe sequential 
changes of the Si dimer and Mn trimer configuration. 
Our results are important for the successful integration 
of Mn atoms into future silicon spintronics devices. 

The filled-state STM image in Fig.Qb) shows Mn 
wires near a step edge of the Si(001) surface. Here, about 
a tenth of a monolayer of Mn was deposited at RT with 
a rate of 55 pm/min. The wires, with typical lengths 
ranging from 5 nm to 50 nm, frequently nucleate at de- 
fect sites or step edges, where clustering of Mn wires is 
observed (dashed ellipse). This reflects the high mobil- 
ity of Mn on the Si(001) surface at RT [21 . Even at 




FIG. 1. (a) Top and side view of the trimer wire structure 
(III) with an atom in the hollow site centered between two 
dimers of a dimer row. Mn atom positions for a dimer wire 
(II) and a wire with one Mn atom at the cave site (I) (we use 
the same terminology as in Ref. [22 ). (b) STM image taken 
at It = 0.15 nA and V s = —2.0 V. Ellipse marks accumulation 
of Mn wires (bright lines) at Si step edge. Circles denote Mn 
clusters. 



these low-temperature growth conditions, a number of 
larger clusters are found (dashed circles). In the follow- 
ing, we take a closer look at two wires, denoted A and B 
in Fig.[TJb) and shown in higher resolution in Fig.[2|a) 
and (d). Interestingly, the two wires have a different 
appearance even though the reconstruction of the Si sur- 
face looks similar in the two images. Two things may 
be noted: First, the Si dimer buckling is pinned near 
the Mn wire. It has been argued [25] that strain along 
the dimer rows can immobilize the Si dimers, which oth- 
erwise move quickly back and forth giving the averaged 
appearance of the well-known Si 2x1 reconstruction at 
RT. From the decay of the apparent dimer buckling (not 
shown) , we estimate that the strain originates at the wire 
and decays over a distance of about 5 nm. The second ob- 
servation is that the Mn wires in filled state images have 
an asymmetric shape with respect to the wire axis with a 
sawtooth-like appearance. In addition to this asymmetry, 
Mn subunits between two Si dimer rows are also tilted. In 
Fig. |2|a) and (d) we highlight this tilt by an open ellipse 
in which the filled circle marks the side with the toothed 
edge. The phase of the Si dimer buckling stays continu- 
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FIG. 2. STM images of two wires A (a) + (b) and B (d) + (e). 
Tunneling parameters were It — 0.2 nA and V s — — 1.7 V 
(a) + (d), —0.75 V (b) + (e). The four wire types that occur in 
our experiments are shown schematically in (c). 



oils across the Mn wire, as indicated by the dashed lines 
in Fig. [2] The appearance of the wires does not depend 
on the wire length, but can sometimes change abruptly 
because of nearby defects. Experimentally, we find four 
types of wires as indicated in Fig.[2jc) with equal prob- 
ability. We reference these in the text using the symbols 
0,0,0 and . The tilt of the wire subunits is al- 
ways the same for a given Si dimer buckling phase (see 
dashed lines). For a given tilt, the sawtooth appearance 
is found to have two realizations with a rotation by 180° 
connecting the two. 

The two panels (b) and (e) in Fig. [2] show the asymme- 
try of the two wires when imaged at small empty state 
bias voltage V s = +0.75 V. The Si dimer atoms closest to 
the Mn wire are again hidden by the wire contrast. How- 
ever, the two next-nearest Si dimers exhibit an atom-like 
appearance similar to that of normal Si dimers when the 
empty dangling bond states are imaged at more positive 
V s . The asymmetry shows up as a dark gap on the side 
of the wire that was smooth in the filled-state images. 
The tilt of the wire subunits is not clearly visible and the 
highest apparent point (marked by x) has moved away 
from the toothed edge towards the cave site. 

Comparing this behavior to that of other metal wires 
on Si [9j [TTJ [13] or to initial stages of Si and Ge growth 
on Si(001) [26], it is clear that Mn is quite unique in 
its asymmetric appearance. All other metal wires are 
attributed to the parallel dimer wire model [see II in 
Fig.JTJa)] and are usually explained by a formation mech- 
anism first described for Al [19]. For these wires, filled- 
state STM images taken under similar conditions look 
completely symmetric with respect to the wire axis, ex- 
cept that in some cases a buckling of the metal dimers 
along the wire is observed [27]. As far as we know, none 
of these metal wires induces significant pinning of the 
Si dimer buckling close to the wire when imaged at RT. 
This suggests that the Mn wire structure creates more 



strain along the Si dimer rows. DFT calculations have 
shown that the hollow site is the most stable adsorption 
site for a single Mn atom on the Si(001) surface [23| [28] . 
For one, two or three Mn adatoms, it is favorable to have 
an atom in the hollow site before the parallel dimer sites 
near the cave site are occupied. This was suggested to 
result in a trimer wire structure as shown in Fig.[TJa) III 
[23] . The additional atom in the hollow site could explain 
the experimentally observed pinning of the Si dimers. It 
is more difficult to explain the observed asymmetric ap- 
pearance of the Mn wire because all three proposed Mn 
configurations in Fig.[TJa) are symmetric with respect to 
the wire axis. However, together with the pinned buck- 
ling of the Si dimers, the trimer wire structure is the only 
one that clearly breaks the twofold rotational symmetry 
present in symmetric metal wires. A rotation by 180° 
keeps the dimer buckling on the surface the same, but 
the configuration of the two buckled Mn atoms near the 
cave site is mirrored. We therefore focus further discus- 
sion on comparing simulated STM images of the trimer 
wire structure at various tip-sample biases with our ex- 
perimental results to identify the source of the unique 
asymmetric appearance of the Mn wires. 

To this end, we carried out state-of-the-art spin- 
polarized DFT calculations [29[ [30° employing the highly 
optimized CPMD code [31]. The PBE ( Per dew-Bur ke- 
Ernzerhof) exchange-correlation functional [32] was ap- 
plied and ab-initio norm-conserving pseudopotentials [33] 
were used. The semilocal pseudopotentials were fur- 
ther transformed into fully separable Kleinman-Bylander 
pseudopotentials [34] , with the d-potential chosen as the 
local potential for the Mn atoms, p for Si atoms, and s for 
H atoms. The electron density was calculated from a 2 x 1 
/c-point grid. The calculated wavefunctions from the T- 
point were analytically extended in z-direction away from 
the surface and integrated over the bias window. Finally, 
the z-value of constant density surfaces was plotted on 
the same grayscale as in the experiment. The results 
are shown in Fig. [3] where the top row shows experimen- 
tal const ant- current images. Here, the sample bias volt- 
age V s is varied from filled- to empty-state imaging from 
left (a) to right (f). The bottom row shows the corre- 
sponding simulated images taking into account a slightly 
smaller effective band gap in the calculations. The im- 
ages show very good agreement for the Si surface where 
the change in contrast for the Si dimers in all images 
is accurately reproduced up to the 2nd Si dimer away 
from the Mn wire. The observed change in appearance 
of the 2nd Si dimers in the empty-state images at low 
bias [Fig.[3jd),(e),(j),(k)] is related to a charge-transfer 
effect from the wire to the 2nd Si dimers. The atom po- 
sitions of these Si dimers [see Fig.|3jm)] exhibit the usual 
buckled behavior with a buckling angle of 18°. This is 
not significantly different from the 18.7° found for the 4th 
dimer pair, consistent with the bare Si value of 19 ± 1° 
in the literature [26] [35]. Nevertheless, the appearance 
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FIG. 3. (a)-(f) Drift corrected STM images as a function of bias (a) V s = -1.75 V, (b) -1.00 V, (c) -0.65 V, (d) +0.65 V, (e) 
+ 1.00 V and (f) +1.75 V. The tunneling current was I t = 74 pA for (c)+(d) and 154 pA for all other images. The calculated 
positions of the Si dimer atoms are marked by the black overlay, (g)-(l) Simulated STM images from DFT (g) Eb — —1.50 eV, 
(h) -0.75 eV, (i) -0.40 eV, (j) +0.50 eV, (k) +0.85 eV and (1) +1.60 eV. (m) Relaxed atomic positions of the 8x4 supercell. 
Arrows to the right of panel (f) indicate the positions of the 1st to 4th Si dimers away from the Mn wire. 



of the 2nd dimers in Fig.[3[d) and (j) is roughly that of 
the 3rd and 4th dimers in Fig.[3|e) and (k) taken at a 
bias which is 0.35 eV higher. This indicates higher filling 
of the Si dimer states near the wire. A similar shift in 
contrast was observed for Ge [26] and Al [10] wires, even 
though in these cases the hollow site is not occupied by 
a wire atom. 

If we look at the wire appearance itself, the most ob- 
vious discrepancy between experiment and simulation is 
found in Fig. (3^d) and (j). The experimental image shows 
a horizontally elongated contrast centered on the hollow 
site, whereas the simulation displays a double-lobed ver- 
tical contrast centered over the highest Mn atom. In 
addition, the filled-state images in the simulation show 
nearly no bias dependance, whereas the experimental im- 
ages exhibit a shift of the highest apparent point in the 
wire subunits indicated by x in Figs.[3|a), (b) and (c). 
We might expect such deviations due to several factors. 
For one, the simulation is done at T = K, whereas 
all experimental images are taken at RT. Additionally, 
our assumption of a flat density of states of the STM tip 
reduces the variations between simulated filled-state im- 
ages, especially at large bias. The periodic boundary con- 
ditions in the simulation further limit the appearance of 
wire asymmetries because the supercell has to be dipole 
free and charge transfer from one Mn wire edge to the 
other may be suppressed. Moreover, even if pure DFT 
simulations are accurate enough to describe the inter- 
action and atomic structure of Mn on Si surfaces [22]- 
l24l [28] , the contraction of partially occupied d-shells is 
often underestimated, resulting in a modified orbital hy- 
bridization and therefore distorted simulated STM im- 
ages. Nevertheless, we show that a closer look at the 



relaxed atom positions allows us to identify the origin of 
the observed asymmetry. 

Despite significant differences between the empty-state 
images at low bias in Fig.[3jd) and (j), the simulation also 
reveals a slightly darker appearance of the 2nd dimers 
above the wire. This arises from a combination of two 
asymmetries, one due to the buckling of the Mn atoms 
along the wire, the other linked to the pinned Si dimers 
near the wire. The situation is shown schematically in 
Fig.JiJa). The two Si atoms (3,4) close to the highest 
Mn atom (B) are found to be 7% further apart than Si 
atoms 1 and 2. Bond lengths 1-4 and 2-3 (see dashed 
lines) are the same on both sides and, with 3.1 A , about 
30% longer than that of normal Si dimers. The bond 
buckling angle is 3.4° for bond 1-4 and 1.8° for bond 2-3 
such that Si atoms 3 and 4 are slightly higher than 1 and 
2. Distance 2-7 is 4 A and identical to 3-8. In contrast to 
this, 4-6 is 3% smaller and 5-1 is 3% larger. This differ- 
ence in proximity explains the higher contrast between 
the two atoms of the 2nd dimer (5,6) above the wire and 
the lower contrast for the 2nd dimer atoms below the 
wire (7,8) in the simulation. Furthermore, it leads to an 
overall darker appearance of the 2nd dimers above the 
wire, in agreement with the experiment. 

The mechanism for the asymmetric appearance relies 
on the relative orientation of the pinned Si dimer buck- 
ling and the buckling direction of the two Mn atoms near 
the cave site. Following our previous notation, the wire in 
Fig. [^a) is a© wire and, together with , has the high- 
est Mn atom to the left of the cave site. Similarly, Q and 
have the highest Mn atom to the right of the cave site. 
Fig.|4jb) and (c) show an example to highlight how the 
buckling directions of Si and Mn determine the appear- 



4 




O Si Mn ■ ' 

i , , > /J 



FIG. 4. (a) Schematic of relaxed atom positions. Circle size 
indicates height of the atoms. The distortion of the Si atoms 
near the hollow site is exaggerated to illustrate the origin 
of the asymmetry, (b) and (c) Successive STM images of 
wire B before full relaxation of the Si dimer buckling. STM 
parameters were I t — 0.2 nA, V s — —2.0 V (b) and V s 
1.7 V (c). Black arrows indicate the slow scan direction. 



ance and atomic structure of the Mn wire as evidenced 
in sequential STM scans of wire B done shortly after Mn 
deposition. A rearrangement of Mn atoms during these 
scans changed the local environment sufficiently so that 
both the Si dimer buckling phase and the buckling direc- 
tion of the Mn atoms along the wire was modified. Start- 
ing at the bottom in Fig.Um), the wire initially has a 
appearance. In the top part, imaging shows instability 
followed by a switch toO for the last two wire subunits. 
This arises from a change in the buckling phase of the 
surrounding Si dimers. Scanning back down in Fig[4^c), 
the Si dimer buckling phase is now continuous over al- 
most the entire image, but the apparent wire asymmetry 
still changes. This points to a switch of the Mn dimer 
buckling between the top O and the middle © part of 
the wire. Both the Si buckling phase and the Mn atom 
buckling change again at the bottom end of the wire, 
where two subunits in the configuration remain. A 
further rearrangement of Mn atoms at the lower end of 
the wire finally converted almost the entire wire to © 
, as presented in Fig.[2jd). The two topmost subunits 
of the wire remained pinned in© and made scanning in 
that area unstable (not shown). This can be explained 
by a pinning of the Mn trimer orientation in opposite di- 
rections at the two ends of the wire, leaving an unstable 
boundary that is not present in other wires even though 
the Si dimer buckling phase is the same along the entire 
wire. 

In conclusion, the unique asymmetric appearance of 
Mn atom wires on the Si (001) agrees best with an ex- 
tended trimer model. The asymmetry arises from the 
interplay between the orientation of the Si dimer buck- 



ling and the buckling of the Mn atoms near the cave 
site. Although many features seen experimentally could 
be reproduced in the simulated STM images, more re- 
fined models are needed to reproduce the exact appear- 
ance of the observed asymmetry. This is of particular 
interest in light of recent predictions of exciting spin and 
transport properties of Mn wires [24], which, if proven to 
be correct, may open the door to develop novel atomic- 
scale spintronic devices in silicon. 
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